ENERGY TRANSPORT IN STOCHASTICALLY PERTURBED 

LATTICE DYNAMICS 
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Abstract. We consider lattice dynamics with a small stochastic perturbation 
of order e and prove that for a space-time scale of order e _1 the local spectral 
density (Wigner function) evolves according to a linear transport equation 
describing inelastic collisions. For an energy and momentum conserving chain 
the transport equation predicts a slow decay, as 1/ y/i, for the energy current 
correlation in equilibrium. This is in agreement with previous studies using a 
different method. 



1. Introduction 

Over the recent years there have been strong efforts to understand energy 
transport in anharmonic chains. A prototypical example is the Fermi-Pasta- 
Ulam chain with pinning potential and a quartic nonlinearity. Its hamiltonian 
reads 

N N-1 

h (9,p) = E + X??} + E { W^+i - <n? + - ?i) 4 } • (!) 

j=-N j=-N 

Here qj is the displacement of the j-th particle from its equilibrium position j, pj is 
the canonically conjugate momentum, m is the mass of the particle. The potential 
has two parts. The term uj^qj confines the j-particle and is therefore referred to 
as pinning potential. The second part depends only on displacement differences. 
It consists of a harmonic nearest neighbor term of strength a% > 0, and the 
anharmonicity P(qj+i — qj) 4 ,P > 0. Energy transport can be studied by coupling 
the end particles, g_jv and q^, to thermal reservoirs at different temperatures. 
An alternative method would be to take first N — > oo in thermal equilibrium 
and to monitor the spreading of the energy when initially deposited close to the 
origin. If f3 — 0, the energy spreading is ballistic and the steady state energy flux, 
jjv, is independent of N. For (3 > and for non-zero pinning one finds, mostly 
based on numerical simulations, that the current ]n = 0(1/ N). This is called 
regular transport. Generically the spreading of energy is then diffusive. On the 
other hand, for u>o = 0, the energy transport is anomalous, ]n = N~ a with a in 
the range from 1/3 to 2/5 (see [11] for a review on the subject). 
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The mathematical analysis of the energy transport in the FPU-chain is a dif- 
ficult task, in particular since non-zero temperature is required. In |4] one of 
the author's proposed to replace the nonlinearity in (CQ) by a stochastic exchange 
of momentum between neighboring sites such that the local energy is conserved, 
with the extension to local momentum conservation being worked out in [21 [3]. We 
refer to section [2] below for a precise definition. The goal of our paper is to under- 
stand in this model the mechanism of energy transport for small noise strength, 
in particular to determine conditions for regular, resp. anomalous, transport. 

Energy transport refers to large space-time scales. In case (3 = the equations 
of motion for (PJ reduce to a discrete wave equation, and semiclassical analysis 
provides a convenient tool. For the harmonic chain at non-zero temperature such 
a program was carried out by Dobrushin et al. in a series of papers [5JEIE!- The 
local spectral density W, in physics parlance the Wigner function, is governed by 
the linear transport equation 

-W(r, k, t) + —u;'(k)—W(r, k, t) = . (2) 



W depends on position reR along the chain, the wave number k £ [— |, |], and 



time t. In our model the dispersion relation u is computed from 

oj(k) 2 = ul + ai (1 - cos(2vr/c)) . (3) 
The local energy density at time t is defined through 

rl/2 

e(r, t) = / W(r,k,t)dk (4) 

J-l/2 

in our units. Thereby ([2]) provides a quantitative description for the energy 
transport in the harmonic chain. 

In [S] the analysis of Dobrushin et al. is extended to a wider class of harmonic 
lattices, including higher dimensions. Mielke [14J covers at great generality the 
case of initial conditions of finite energy, see also [9]. 

The starting point of our investigation is the observation that the local spectral 
density should also be the appropriate quantity when stochastically perturbing 
the lattice dynamics, provided the strength of the stochastic term is appropriately 
adjusted to the semiclassical limit. This program indeed works out and, as our 
main result, we will prove that the stochastic exchange gives rise to a linear 
collision term in (j2J) and the transport equation is modified to 



^-W(r, k, t) + —uj'Vk)^-W{r, k,t) = [ R(k, k') (W(r, k\ t) - W(r, k, t)) dk . 
at 2w or J 

(5) 

with an explicit transition kernel R(k, k'). 

We will give in detail the proofs in the one-dimensional case with a stochastic 
dynamics that conserves total momentum and total energy (as in [2]). These 
proofs can be straightfowardly extended to lattices of dimension d > 2, and to 
stochastic perturbations that conserves only energy (as in [4]). 
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The limit transport equation is easily analysed and is a convenient tool to study 
energy transport. For the one dimensional model with conservation of energy and 
momentum, one finds that for uq 7^ the energy transport is diffusive while for 
no pinning, uo = 0, the spreading is superdiffusive, in accord with the results in 
[21 [3]. The super diffusion process is governed by a Levy process of index | (see 
|10j). If ujq > one finds regular energy transport, the same holds for dimension 
d > 3 and for noise conserving only energy, in accordance with a direct approach 
using the Green-Kubo formula [21 [31 H] . 

Our paper is organized as follows. In sections [2] to H] we derive equation (j2J) 
including the collision term. In section [5] we explain how our proof generalizes 
to higher dimensional lattices. We then discuss the energy transport in detail 
(section [6]). 

2. The model (One-dimensional case) 

To develop the necessary techniques, we consider first the case of a one-dimensi- 
onal chain, as in ([T]) with (3 = 0. In view of higher dimension it is convenient 
to label the particles by y G Z. The chain is infinite and we allow for a general 
harmonic coupling. The phase space is (R x IR) Z and a configuration at time t is 
denoted by {q y (t),Py(t)} ye z- The Hamiltonian of the system is given by 

h(p, q) = \^p 2 y + \ Yl a ( y ~ y')wv'i ( 6 ) 

where we use units such that m = 1. We denote with v(k), k e T = [0, 1], the 
Fourier transform of a function v on Z, 

v(k) =J2e~ 2mkz v(z), (7) 

and with f(z), z G Z, the inverse Fourier transform of a function / on T, 

/(*)= / dke 2 ™ kz f(k). (8) 

The function u>(k) = \J a{k) is called dispersion relation. 
We assume a(-) to satisfy the following properties: 

Assumption 1. 

• (al) a(y) 7^ for some y 7^ 0. 

• (a2) a(y) = a(—y) for all y e Z. 

• (a3) There are constants C±, C2 > such that for all y 

Hy)\ < C x e~ c ^. 

• (a4) We require either 

— (pinning): a(k) > for all k G T 

- (no pinning): a(k) > Wk ^ 0, d(0) = 0, d"(0) > 0. 
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Assumptions (a2), (a3) ensure that a is a real analytic function on T. Note 
that the potential depends only on the displacements differences iff u)(0) = 0, 
hence the two cases in (a4). In the pinned case 00 is strictly positive and thus 
real analytic on T. In the unpinned case, the condition (a4) says that u(k) = c\k\ 
with c > for small k, to say uj is a regular acoustic dispersion relation. 

We consider the Hamiltonian dynamics weakly perturbed by a stochastic noise 
acting only on momenta and locally preserving momentum and kinetic energy. 
The generator of the dynamics is 

L = A + e 1 S (9) 
with e > 0, where A is the usual Hamiltonian vector field 

A = J2 p y d iv ~ J2 a ( y ~~ y')% ,d vyi (10) 

y& y,y'eZ 

while S is the generator of the stochastic perturbation. The operator S acts only 
on the momenta {p y } and generates a diffusion on the surface of constant kinetic 
energy and constant momentum. S is defined as 

where 

Yz = (p z ~ + (Pz+i - Pz-i)d Pz + (p z -i -p z )d Pz+1 

which is a vector field tangent to the surface of constant kinetic energy and of con- 
stant momentum for three neighbouring particles. As a consequence energy and 
momentum are locally conserved which, of course, implies also the conservation 
of total momentum and total energy of the system, 



£^ = 0, SH = 0. 



The evolution of {p(t),q(t)} is given by the following stochastic differential 
equations 

alq y = p y alt, 

dp y = - (a* q) y alt + —A(Ap y + p y „i + p y+ i)dt 

6 (12) 



+ V y Yl ( Y v+kPy) dw v+k W • 

' fe=-l,0,l 



Here {w y (t)} y£ z are independent standard Wiener processes and A is the discrete 
laplacian on Z, 

Af(z) = f(z+l) + f(z-l)-2f(z). 

To study the local spectral density it is convenient to introduce the complex 
valued field ip : Z — > C defined as 

^(y,t) = ^={(u*q) y (t)+tp y (t)). (13) 
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Observe that \ip(y) | 2 = + | J2y'ez a (,y~y') ( ly<ly' = e y * s ^ ne energy of particle y 
and conservation of total energy is equivalent to the conservation of the £ 2 - n orm. 
For every t > the evolution of if> is given by the stochastic differential equations, 



dip(y,t) = - i(u>*if>)(y,t)dt+-ej/3 * (ip - ip*)(y,t)dt 

+ \M E (Y y+k \{ii)-r){y,t))dw y+k {t), 

" fc=-l,0,l 



(14) 



where (3 is defined through 

09 * /)(*) = U(4/(z) + /(* - 1) + f(z + 1)). (15) 

D 

3. WlGNER DISTRIBUTION AND THE BOLTZMANN PHONON EQUATION 

Given a complex valued function J on K x T, we define onKxZ 

J{x,z) = [ dke 2mkz J(x,k). (16) 

We also define oeRxT 

J(p,k)= [ dxe~ 27ripx J(x,k). (17) 
Jr 

We choose a class of test-functions J on 1 x T such that J(-,k) G «S(R, C) for 
any fceT. 

Let us fix e > 0. We denote by ( • ) the expectation value with respect 
to a family of probability measures on phase space which satisfies the following 
properties: 

(bl) {ifj{y)) e = 0, \/y G Z; 

(b2) (WMv)) e = 0, Vy,|/GZ; 

(b3) sup £>0 £ (|M| 2 2 ) e < # for some AT > 0. 

Observe that, since (||'?/'|| 2 ) e = (H) is the expectation value of the energy, we 
are considering states with an energy of order e^ 1 . 

For every if> we define the associated Wigner function in the standard way, see 
[TJ [121 H31 US] , and integrate it against the test function J. The bilinear expression 
in if) is averaged over < • > £ . For simplicity we call the averaged Wigner function 
simply Wigner distribution and denote it by W £ with e the small semiclassical 
parameter. Thereby we arrive at the following definition 

( J > W£ ) =\ E Wmv)) E J dke 2 ^y'-y\j(e(y' + y)/2,ky 

=\ E Wrm)J«y' + y)/2,y-yr, 

where J G 5(1 x T). By condition (b3) this distribution is well defined, see 
proposition [13] in the appendix. 
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Proposition 2. Under the assumption (b3) for every test function J G <S(R x T), 
there exist constants Ki,K 2 such that 

\(J,W 6 )\ <Ki [ dp sup|J(p,Jfe)| < oo, 
|(J,W £ )| < if 2 ^ sup|J(a;, z)\ < oo 
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/or every e > 0. 



Remark 3. Notice that W £ is well defined on a wider class of test functions than 
S(M. d x T d ). In particular we can take J(x,k) = J{k), a bounded real valued 
function on T, and by ( flgj) we have 

(j,w £ ) = £ -^dk(mk)f) e j(k), 

while choosing J(x,k) = J(x), a bounded real valued function on R, we have 



Let us start our dynamics with an initial measure satisfying conditions (bl), 
(b2), (b3). We want to study the evolution of the Wigner distribution W e on the 
time scale i.e. we define for J G «S(R x T), 

= \ E W> <A) V(y, t/e)) E f dk e 2 ^y'-^J(e(y' + y)/2, k)* . ^ 

Observe that since the dynamics preserves the total energy, the condition e (H^ll) 
< K holds at any time and, by proposition [T3J the Wigner distribution is well 
defined at any time. 

According to remark [21 if we choose test functions depending only on k, then 
we obtain the distribution of energy in fc-space. It turns out that in the limit as 
e — > 0, this distribution converges to the solution of the homogeneous Boltzmann 
equation, and this will be our first result. We define the distribution S £ (t) on T 
by 

(j,S%t)) = (J,W%t)) = \f dk (mk,t/e)\ 2 ) £ J(k) (21) 

for any bounded real valued function J, and introduce the collision operator C 
acting on «S(R x T), 



CJ(x, k) = / dk'R(k, k') (J(x, k') - J(x, k)) . (22) 
Jt 

with 

4 

R(k, k') = - (2 sin 2 (27rA;) sin 2 (vrA;')+2 sin 2 (27rfc') sin 2 (vr£;) 

— sin (2nk) sm 2 (2nk')). 
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Theorem 4. Assume that the initial measure satisfies conditions (bl),(b2), (b3), 
and furthermore S £ (0) converges to a positive measure So(dk) on T. Then £ e (t) 
converges to £(t,dk), the solution of 

d t {J,£(t))= 1 (CJ,£(t)) (24) 

for every bounded function J : T — > M. 

In order to prove the next theorem, the full inhomogeneous equation, we need 
an additional condition on the initial distribution in the unpinned case (a(0) = 0): 

(b4) In the case of no pinning we require 

lim EE - / dk (\ip(k)\ 2 ) = 0. 

B-o ~o 2 y, fe | <fl 

Condition (b4) ensures that there is no initial concentration of energy at wave 
number k — 0. This condition can be omitted for a dispersion relation u which 
is analytic on T (as in the pinned case). 

Theorem 5. Let Assumptions (bl-b4) hold and assume that W e (0) converges 
to a positive measure fi (dx,dk). Then, for all t e [0, T], W £ (t) converges to a 
positive measure fi(t) = dx, dk), which is the unique solution of the Boltzmann 
equation 

d t {j^{t)) = ^{J(k)d x J 1 n(t))+ 1 {CJ^{t)) (25) 
with initial condition ^(0,dx,dk) = fi (dx,dk). 

In (1251) (J, fi(t)) denotes the linear functional L T J(x, k)*n(t, dx, dk). 

Observe that the kernel R of (|23l) is non-negative, symmetric, and is equal to 
zero only if k = or k! = 0. Moreover, it is easy to see that 

/ dk' R(k,k') = -p(k), (26) 
./ 

where /3(k) is the Fourier transform of the function (3 defined in (1151) . Thus 
the Boltzmann equation (1251) can be interpreted as the forward equation of a 
Markov process (X(t),K(t)) on R x T for the dynamics of a particle, which in 
the context of lattice dynamics is called phonon. The phonon with momentum k 
travels with velocity u'(k)/27r and suffers random collisions. More precisely K(t) 
is an autonomous reversible jump Markov process with jump rate R, while the 
position X(t) is determined through 

f l 1 

X(t) = X(0) + / —u'(K(s)) ds. 
Jo 2tt 
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4. Proof of theorems EE and [51 

4.1. Relative Compactness of the Wigner distribution. Existence of the 
limit of the Wigner distributions W e (t) will be established as in [T]. The limit 
distribution W(t) is non-negative, as it is proved in [T2], |T3], i.e. for every 
J E S(R x T, C) one has 




Let us introduce the space A of functions J on M d x T d such that 

|| J\\a = 2_] SU P I J{ x i z )\ < °°> (28) 

where J is defined in (|T6j) . The following lemma shows that if the distributions 
W £ are uniformly bounded in A', the dual space to A, then at every time t one 
can choose a sequence £j — > such that W £l converge in the *-weak topology in 
A' to a limit distribution W(t). 

Lemma 6. There exists a constant C > independent of t such that We > 

\\W £ {t)\\ A ,<C. 

Proof. For every J e *4 

<4W 8 (t)> = f (*M(v))J{e{v + \f)l%V-lfY. 

Then, using Schwarz inequality and Assumption (b3), 

\(J,W°(t)) e \ < |^<|^)| 2 )JJ|U < K\\J\\ A . 

□ 

4.2. Proof of theorem [41 We consider a class of test functions J depending 
only on k G T. In particular we choose J real valued and bounded. Recall the 
definition 

(j,£ £ ) = yjk{\m\\j(k), 

which is well defined since |(J, £ e )| < |i^sup fcgT |J(A;)|. 

The evolution of the distribution S £ (t) is determined by Ito's formula, namely 

d t (J,£%t)) = - I dke-\L\^Ktle)^) e J{k) , 

where 

L\m\ 2 = A\m\ 2 +e 1 s\m\ 2 

and A, S are respectively defined in ( flOl) . ( fTTI) . We have 

A|^(A;)| 2 = [A$(k)*]$(k)+$(ky[A$(k)], 
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where by direct computation 

Aif}(k) = £ e~ 2mky A^j{y) = -i^2 e~ 2mky {^ * if>) (y) = -iuj{k)jj{k) (29) 

and thus ^^(/c) | 2 = 0. For the stochastic part, since S is a second order operator, 
we have 

s\m\ 2 = [sj(kyw(k) + j,(k)*[sj(k)\ + 

where by direct computation 

y&L yeZ (30) 

with (3 defined in ( 115ft . Thus 

(s$(ky)ii>(k)+$(k)*(sj>(k)) 

where 

4 

(3(h) = -- sin 2 (7rA;)(l + 2 cos 2 (vrA;)). (31) 
o 

Finally we have to compute \ Ylzez\^^{k)*]\X z '^'{k)]- It holds 

J2[YMky][YMk)] = £ e 2 ^'^ £[r^(y')1[^(2/)], (32 ) 

where 

= K^tKiOltWCi/)] + [Yy^yTWyHy)] + [V^l^KV^fe)]- 

This expression is explicitly computed in the appendix, see eq. (IHTj) . By inserting 
it in (1321 we get 

£ e 2 ^^-^£[r^(y01[n^(y)] 

= cos(4ttA;) £(2p^ + i - p y p y+2 - p 2 y ) 

y&Z 

+ cos(2ttA;) ^2(2p y p y+ 2 - 2p 2 y ) + £(-2^^+! - p y p y+2 + 3p 2 ), 

y ez y ez 



10 GIADA BASILE, STEFANO OLLA, AND HERBERT SPOHN 

which is equal to 

/ d£ |p(0| 2 (cos(4vrA;)[2cos(27rO - cos(4<) - 1] 
Jt 

+ 2cos(27rA;)[cos(47r£) - 1] + [3 - 2cos(2tt£) - cos(47r£)])- 
Finally, after some trigonometric identities and using the relation 

\p{k)\ 2 = \{\m\ 2 + m-k)\ 2 - m^(-k) - twfc-k)*), 

we get 

where R(k,£) is given by ( 1231) . 

Since f T d^ R(k,£) = —/3(k), we can write 

simi 2 = c\m\ 2 - lc$(k)#{-k) + #(ky${-ky), 

where C is the operator defined in (122]) . i.e. 

Cf(k)= [ d£R(k,0{f(0-f(k)). 
Jt 

The evolution of £ £ {t) is given by 

d t (J, S £ (t)) =-Y^Jdk (mk, t/s)\ 2 )CJ(k) 

-^l dk \[{m)^-k)){t/e)) £ + ((m^(-kY)(t/e)) £ ](CJ)(k). 
Defining the distribution Y £ (t) on T through 

<^ £ W) =\ E W,t/e)Hy,t/e)) E f dke 2 ™ k ^J(k) 

y,y'eZ Jj 

= -J^dk ([m^(-k))(t/e)) e J(k), 

we can rewrite the evolution equation as 

d t (J,S e (t)) = 1 (CJ,£ £ (t)) - ^((CJ,Y £ (t)) + (CJ,Y £ (t)*)). (33) 

This is not a closed equation for £ £ (t). However we expect that in the limit e — > 
the terms containing the distributions Y £ (t), Y £ (t)* disappear. In order to prove 
it, let us consider the evolution of the distribution Y £ (t) on the kinetic time scale. 
Calculations are similar to the previous ones, but with the difference that now 
A[ip(k)ifj(—k)} 0, and precisely 

A${k)${-k)\ = -2iu{k)${k)${-k). 
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We arrive at 

d t (J, Y £ {t)) = J, Y £ {t)) + 10 J, Y £ {t)) 

+ l((cj,Y%t)) + (c.j,Y%ty)) - lifavityy-rtwit)). 

Observe that by integrating eq. (I34|) in time, we obtain 



(34) 



lim 

£^0 



dt (uJ,Y £ (t)) 







for every bounded function J. In particular, since by item (i) of lemma [T 

R(k, k' 



sup 

fcGT u[k) 



< OO, 



we can choose a function uo CJ with J bounded and obtain 



lim 

£^0 



dt (CJ,Y £ (t)) 



0. 



In the same way we have lim e ^o | J* * dt (CJ, Y £ (t)*^ | = and any limit distribu- 
tion £(t) of £ £ {t) solves the equation 



(j,S(t)) = (J,£(0)) + 1 / ds(CJ,£(s)) 

Jo 

for every bounded real valued function J. □ 

4.3. Proof of theorem [51 Now we will give the proof of (|25|) for test functions 
J G 5(lx T, C). The main difference to the previous case is that the Hamiltonian 
part of the generator contributes to the evolution of W £ (t), resulting in a ballistic 
transport term. In order to control this term, we need to ensure that there is no 
mass concentration at k = for every macroscopic time t G [0, T] with T > 0. 
This is stated in the following lemma. 

Lemma 7. Let assumption (b4) hold. Then for every t G [0,T] 



lim lim 

p^o £^o 2 



k\<p 



dk(mk,t/e)\ 2 ) t 



0. 



Proof. We use the evolution equation f[33l for J p (k) = l[_ pp ](A;). Since 
\CJ p (k)\ < Cl (2p + J p (k)) and (\CJ p (k)\,\Y £ (t) + Y £ (tY\) < (\CJ p (k)\,£ £ (t)),we 
obtain the bound 

(J P ,S £ (t)) < (J p ,£ £ (0))+c 2l [ ds (\CJ p \,£ £ (s)) 



< (J P ,£ £ (0)) + c 3 j(2 P Kt + J ds (J P ,£ £ (s)) 

where K is the bound on the total energy from condition (b3). Then by Gronwall's 
inequality 

(J p ,£ £ (t)) < (2pK + (j p ,£ £ (0)))e c ^, 
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where, by assumption (b4), \im £ ^ (j p , £ e (0)} — ► for p — ► 0. □ 

4.3.1. Proof of theorem^ For every J G S(R x T, C) the evolution of the distri- 
bution W £ (t) on the kinetic time-scale is given by 

d t (J,W s (t)) 

= (e/2) d t ^(y',t/e)^(y,t/e)) £ f dk e^'^J{e{y + y')/2, k)* 
= {e/2) e- l {L[i,{y'Y^{y)]) £ f dk e 2 ^'^ J{e{y + y')/2, k)\ 

where 

= A^{y'y^y)}+e 1 S[^y'r^y)] 

and A, S are defined in (fTOj) . (TTTj) . respectively. We start by computing the 
evolution determined by A, the Hamiltonian part of the generator. Using the 
representation of the Wigner distribution in Fourier space we get 

e 

2 



T 



dp I dk s^A^ipik - sp/2)*tjj(k + sp/2)) e }J(p,ky 

dp / dk(tp(k - ep/2)*^{k + ep/2)) £ 
x e~ 1 [u(k + ep/2) - cu(k - ep/2)) J(p, k)* . 



e 

—i— 
2 



Now we prove that one can replace e l [u(k + ep/2) — u(k — ep/2)} with u'(k)p 
in the last expression. For every < p < 1/2 



dp / dk{4){k-ep/2)*i){k + ep/2)) e 
Jn 

x (e~ x [u{k + ep/2) - u(k - ep/2)] - co'(k)p) J(p, k)* = I* (p) + I* (p), 



where 

H (n) = 



I E >{P) 4 / dp I dk (i>(k - ep/2)^(k + ep/2)) e 

\k\>p 



x (e~ l [oj(k + ep/2) - u(k - ep/2)} - u'{k)p) J{p, k)* 



e 



I<(P) =7j dp dk (^(k - ep/2)^(k + ep/2)) £ 

1 JR J\k\<p 

x (e~ l \uj(k + ep/2) - uj(k - ep/2)} - co'(k)p) J(p, kf . 
Using Schwarz inequality and points (i), (ii) of lemma [H] in the appendix 

\r K {p)\ < [ dp (C\p\ + 11^11^) sup \J(p,k)\(~ [ dk(\^(k)\ 2 ) E ) 



<c £ -i dk(\m\ 2 ) £ 

Z J\k\<P 
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and, by lemma [3, 

limlim \li(p)\ = 0. 
p^Oe^O 1 <v n 

To compute />(p) we split it into two parts, 

I E >{P) =\\ d P I dk (4>(k - ep/2)*i>{k + ep/2)) e 

^IR,e|p|>p J\k\>p 

x (e~ l [u(k + ep/2) - cu(k - ep/2)) - tu'(k)p) J(p, kf 
+ ^ / dpi dk (i>(k-ep/2)*i>(k + ep/2)) £ 

^ Jw.,s\p\<p J\k\>p 

x (e~ l [uj{k + ep/2) - uj(k - ep/2)} - u\k)p) J(p, k)*. 

Again we use Schwarz inequality and points (i), (ii) of lemma fl4l to show that 
the first term on the RHS is negligible, since for all p > it is bounded by 

K dp (c\p\ + || Vc^Hoo) sup | J(p, k)\, 

J\p\>p/e keT 

which tends to as e — > 0. 

For the second term on the RHS we use the point (iii) of lemma [T3J since 
\k\ > p, e\p\ < p implies \k\ > e\p\, and for all p > we get 

e 
2 



dp / dk (i(){k - ep/2)*ip(k + ep/2)) e 

\p\<p/e J\k\>p 



x (e~ l [oj{k + ep/2) - u(k - ep/2)} - u'(k)p) J(p, k)* 

< K I dp e — |p| 2 sup | J(p, k)\, 
Jr P keJ 



which tends to as e — > 0. Then we have 

e 
2 



\ I dp I dk e- x A[{ij){k - ep/2)*ip(k + ep/2)) ] J(p, jfe)" 



r 



| / dp I dk {i,{k - ep/2fi,{k + ep/2)) e {~i p)uo'{k)j\p,k)* + 0{e) (35) 
1 



. Vco>V r J, W £ (t) ) +0(e). 

2lX \ / £ 

We have to compute 

(e/2) (Smv'n(y)]) E [ dke 2 ^y'-^J(e(y + y')/2,ky 

= \ E (smy'my)})J«y + y')/2,y-yT- 

y,y'eZ 

Since S is a second order operator, we have 

sWY^y)} = W)*s^(y) + [S^y'TMy) + 1 YfrMW^to)], 
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where by direct computation 

S1>(y) = ^S Py = * (r - ^){y) (36) 

and 

iKvWCv) + [S^y'TMy) = 4(^OT*p)»-^(y)(^*pV)- 
Integrating by parts and using the symmetry of (3, we can rewrite 

7 (e/2) £ [{4,{ y ySi,{ y )) £ + {{S4,{ y y)4,{y)) £ ]j{e{y + y i)/2, y - y >y 

y,y>& 

i 



^/ 2 ) E wrp 2 /) £ E^) J >^+^+^/ 2 ^-^+^)* 

7 (5/2) £ * ^( y ) JV > e ^ / g( z )j( e ( y + y '_ z )/2,y-y' + z)*. 



Using the energy bound and the properties of the test functions J, one can write 
the first term on right hand side as 

i 



y,y'eZ * zeZ 
= ^/ 2 ) E ^WTPv) e I dke 2 * ik ^P(k)J(e (y + y')/2,ky 



z 

-' 1 

+ 0(e) 

and the same can be done for the other term. Finally we obtain 

(s/2) [Wysm) E + {{s^{yr)m)}J(<y + v')/^y-y'T 

= (e/2) V -4(^(yO*ftf-^(yW> e / d*e^^(A0J(e(v + vO/2,*^ 

+ 0(e). 

(37) 

About the other term in (l36|) first observe that X^zP'W'O/')*] \Yzi>(y)] is just 
a finite sum for any y, y' fixed and, computing it explicitly and identifying terms 
that differ by translations, see details of the computation in appendix I7.ll one 
obtains 



l -{e/2) {[Yz4>{y'yWMy)})J{e{y' + y)/^y-y'y 

y,y',z& 

2 

= (e/ 2 )E E a(z,u)(p y p y+z ) £ J(ey,uY + 0(e), 
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where a(z,u) = a(—z,u) = a(z,—u) = a(u,z) and is given by fl82|) . We can 
rewrite it as 

2 

( £ / 2 )E E a(z,u)( PyPy+z ) £ J(e(y + z/2),uy + 0(e) 

j/gZ z,u=— 2 

= (e/2)J2Yl ot(z,u)( PyPy+z )J(e(y + z/2),uy + 0(e), 
where we put a(z, u) = if \z\ > 2 or \u\ > 2 and by changing variables we obtain 



(e/2) (PyPy')eJ2 a( -y' -V^VW + v)/%u)* +0{e). 

Defining 



(39) 

y,v'' 



R(k, k') = J2J2 e ' 2nikZe ~ 



-2-nik'u , 

we can rewrite ([391 as 



(e/2) ]T < W > £ / dke 2mk ^ I dk'R(k,k')J(e(y' + y)/2,kr + 0(e), 

y,y'eZ J J 

where direct computation gives 
2 

R(k, k') = - (3 - 2 cos(27rfc) - cos(4ttA;) - 2 cos(2vrA;') + 2 cos(2vr(A; / + 2fc)) 
3 

- cos(4ttA;') + 2 cos(27r(2A;' + k)) - cos(2tt(2A;' + 2A;))) 

4 

= - (2 sin 2 (27rfc) sin^TrA/) + 2 sin 2 (27rfc') sin 2 (7r£;) - sin 2 (27rfc) sin 2 (27rfc')) , 
3 

which is the kernel defined in ( 1231) . Using the relation J T dk'R(k, k') = —$(k) and 
(PyPy') £ 

= \[WY4>{y)) £ + (VW(i/)*>J - \lW)i>{y)) e + Wn(yT) £ } 
■j=w(v r )*Pv-4>(y)Pv') e 

= wrm) £ - \[W)m) e + Mvowu 

we obtain 

l\ E (S(«))) £ W + !/)/^'-!/) 

!/.l/'GZ (40) 

= 7 (CJ, - |«CJ, Y £ ) + (CJ, F £ *)) + O(e), 
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where the collision operator C is defined in fl22|) and the distributions Y £ {t), 
Y e (t)* are defined as 



(41) 



<j,F £ (t))=(e/2) Wmy))J«y + y')/2,y-yT, 

y,y'& 

(j,Y%ty) =(e/2) m y rm*)eMy + y')/^y-yr 

y,y>eZ 

for every J G «S(R xT,C). 

The evolution of W £ {t) is not a closed equation, 

d t (j, W e (t)) ={{VuV x J), W £ (t)) + j{(CJ), W £ (t)) 

- ^((CJ),Y £ (t)) - l((CJ),Y%ty) + 0{e). (42) 

However we expect that in the kinetic limit e — > the terms containing the 
distributions Y £ (t),Y e (t)* to disappear. To prove this, we consider the evolution 
of Y £ (t). Again by Ito's formula 

d t (j,Y%t)) = ( e /2) d t (W)^(y))J«y+y')/^y-y'r 

= (e/2) e-\Lty{y')^{y)]) £ J{e{y + y')/2,y- y y 

y,y'eZ 

= (e/2) e- l {A[^y')^y)]) £ J{e{y + y')/2,y~y'y 
+ 7 (e/2) Yl (SMy'My)]) £ j(6(y + y')/2,y-yy. 

y,y'&Z 

For the stochastic part of the generator, by a similar computation as above, we 
obtain 

7(^/2) (SW(y'M y )}) £ J(e(y + y ')/2, y ' -y) 

y,y'eZ 

= l((CJ),Y £ ) + 1{{CJ),Y") + l(0J),Y £ ) - 1{0J),Y**) (43) 

-l[((CJ),W £ ) + ((CJ),W £ *)}+0(e). 

To compute the Hamiltonian contribution to the evolution of Y £ (t), we use the 
representation of Y £ (t) in the Fourier space and get 



{e/2) dp dke~ 1 A[(i)(k-ep/2)i)(k + ep/2)} £ \ 

= -2ie-\e/2) dp dk e _1 (^(A: - ep/2)^{k + ep/2)) 



(oj(k + ep/2) + u(k - ep/2))J(p, k)< 
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where with similar arguments as above one can replace u>(k + sp/2) + ou(k — ep/2) 
by 2uj(k), with an error of order e 2 . Then we arrive at the following equation for 
the evolution of Y £ (t), 



2% 



d t (j, Y%t)) = ~((u,J), Y £ (t)) + 1((CJ), Y £ ) + 1(0J), Y £ ) 
+ 1((CJ),Y**) - |((/3J),Y £ *) -1[((CJ),W £ ) + ((CJ), W £ *)] + 0(s). 



1 



(44) 



After time integration, we obtain, for any J G <S(R x T) 



lim 



dt ((coJ),Y £ (t)) 



0. 



Observe that R(k, k')/uj{k) G C°°(T/{0}) and, using item (i) of lemma fT4l in the 
appendix, 

R(k,k') 



sup 



< oo. 



Then equation (1441) holds for any function 

R(k,q) 



u{ky l CJ(y,k) 



dq 



u(k) 



(J{y,q)-J(y,k)) 



with J G «S(R x T) and consequently 



lim 



dt ((CJ),Y £ (t)) 



0. 



5. Extension to dimensions d > 2 

We consider a particular generalisation of our model to d dimensions, d > 2. 
The perfect lattice is 7L d . Deviations from the equilibrium position y G Z d is 
q y G M. d and p y denotes the corresponding momentum. Thus the phase space is 



x 



The Hamiltonian of the system is given by 

H(p, q) = - p y + o ~ y ') q y ' V- 



(45) 



y.y 



For simplicity the couplings a are taken to be scalar. In general, a would be a 
d x d matrix. We denote 



><■<) = £ 



-27rik-z 



«(*), /(*) 



<ik e 



27rik-z 



/(k). 



(46) 



We assume «(•) to satisfy the following properties: 

Assumption 8. 

• (al) «(y) 7^ /or some y / 0. 

• (a2) a(y) = at(-y) /or a// y G Z d . 
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• (a3) There are constants C 1 , C 2 > such that for all y 

Ky)|<Cie- C2|y| . 

• (a4) 

— (pinning): a > on T d , 

- (no pinning): a(k) > for all k ^ 0, d(0) = ; Hess(a(0)) is 
invertible. 

The dynamics is determined by the generator L = A + ejS with 

A =Y. p y' d *y- E a (y-y')qy'-dp y - (47 ) 

ygz d y,y'ez d 
5" is defined through the vector fields 

x% = U - Pi)(d Pi - d p - (pi - pi)(d pi - d^, 

according to 



S 



^3T)E E Ky+eJ 2 = ^^Ty E E( X ^) 2 > (48) 

yeZ d i,j,k=l y.^ezd, «j'=l 

l|y-*ll=i 



where ei, . . . , is the canonical basis of Z . As in the one-dimensional case 

S Py = > 5iJ = °- 

yei, d 

Note that now it suffices to couple nearest neighbors. 

The evolution of {p(t),q(t)} is given by the following stochastic differential 
equations 

dq y = p y dt, 

dp y = — (a * q) y dt + 2^7 Ap y dt 

d (49) 

+ E E Mi*) 

z6Z d , lj = l 

ll*-y||=i 

for all y G Z d . Here {w^ y = w^ z ; z,y G Z d ; 2, j = f, . . . ,d; ||y — z|| = 1} are 
independent standard Wiener processes. 

As before, we define the complex valued vector field ij> : Z d — > C d by 

^(y> *) = ^ ((^ * q)y(<) + ^Py(O) (50) 



with the inverse relation 



p y (t) = ^=(r-^)(y,t). (51) 
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Observe that |-0(y)| 2 = e y , the local energy ay y. For every t > 0, the evolution 
of i\) is given by the stochastic differential equation, 

di>(y, t) = -i(u* T/>)(y, t)dt + X -e^ * (V - ip*) (y, t)dt 

+ E E(^(^-^)(y,i))^w, (52) 

l|y'-y||=i 

where (3 is determined through (/3 * /)(z) = A/(z). 
Given a function J on M d x T d , we define 



J(x,z) = / dke 42 * z J(x,k) (53) 

on R d x Z 11 . We also define 

J(p,k)= / dxe- &p 'V(x,k). (54) 



We choose a class of test-functions Jonl d x T d such that J G S(R d x T d ,M d ), 
where M rf is the space of complex d x d matrices. 

Fix e > 0. We introduce the complex valued correlation matrices 

<^(yT ® VKy)} £ , <^(y') ® V>(y)> £ , (55) 

where ( • } £ denotes the expectation value with respect to a probability measure 
on phase space which satisfies the following properties: 

(cl) <V(y)> £ = 0, Vy G Z d , 

(c2) (^(y')®^(y)> e = 0, Vy,y'GZ d , 

(c3) <IIV'l| 2 > £ = <E z ^IV'(z)| 2 > £ <^. 

Observe that, since (||V'l| 2 ) e — (-^) £ > we are considering states with an energy of 
order e~ d . We define the matrix-valued Wigner distribution W e as 

(J,W £ ) 

= (e/2) d E <^(y')^(y)> £ / rfke^ k -^'^J,,(,(y' + y)/2,A;)* 

y,y'ez^j=i jTd (56) 

y,y'£Z d i,j=l 

with J G 5(lR d x T d , M. d ). The evolution of the diagonal terms of the distribution 
W £ on time scale e _1 t is determined through 

d 

(j,w%t)) = (e/2) d J2 E <^(yW^(y,^)) £ 

i=l y,y'eZd (57) 

x / rfke^ k -( y '- y )j i (£(y' + y)/2,A;)* 
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for J, G S(M. d x T d ), i = l...,d. 

Observe that since the dynamics preserves the total energy, the condition 
£ d (ll^ll) < K holds at any time and, by proposition [T3l the Wigner distri- 
bution is well defined at any time. On this time scale the diagonal terms of 
the distribution W £ (t) converge in a weak sense to a (vector valued) measure 
/i = {fii(t),i = 1, .., d} on M. d x T d which satisfies the following Boltzmann equa- 
tion. For any vector valued function J G S(M. d x T d , C d ), 

(j,fi{t))-(J,fi{0)) = ^- [ ds ((Vu-V x J,(i(s))+ 1 (CJ,v(s))), (58) 

Z71 Jo 

where (J, denotes the scalar product ^^ =1 / RdxTd Jj(x, k)*/Xj((ix, (ik). The 
collision operator is given by 

(CJ) i (x,k) = ^ T J^' R(KV)(Jj(*,V) - Ji(x,k)), (59) 

where the kernel R : T d x T d — > R has the following expression, 



i?(k, k') = 16 ^ sin 2 (7r^) sin 2 (7r^). (60) 
e=i 

As in the one-dimensional case, in order to prove the inhomogeneous Boltzmann 
equation fl58|) . we need an additional condition on the initial distribution in the 
unpinned case (a(0) = 0), which ensures that there is no initial concentration of 
energy at k = 0: 

(c4) In the case of no pinning we require 



lim Hm( £ /2) d / dk <|V>(k)| 2 ) £ = 0. 

K—0 e^O J\k\<R 



\k\<R 

Now we state the precise theorem. The proof is analogous to the one-dimensional 
case. 

Theorem 9. Let Assumptions (cl-c4) hold and assume that W e (0) converges to 
a positive vector valued measure fiQ. Then, for all t G [0, T], W £ (t) converges 
to a positive (vector-valued) measure p,o(t) which is the unique solution of the 
Boltzmann equation 

d t {J,n{t)) = ^{Vu ■ V X J,M*)> +7<CJ,/i(i)> (61) 
with initial condition Ho(t). 

As in the one-dimensional case, the Boltzmann equation has a probabilistic 
interpretation as the forward equation of a Markov process. We consider the 
Markov process 

(X(f),K(t),i(t)). 
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By (1591) . the jump rate from (i, k) to (j, dk.') is given by 

v Ki {j,dk!) = ^-(1 -8 id )R(k,k')dk', = l,...,d. 

Transitions between states with the same index i are forbidden. The total collision 
rate is 

<Mk) = V/ u k ,(j,dk')= [ dk'i2(k,k'), 

j=i ->^ d ^ d 

i = 1, d, which does not depend on i. Explicitly 

d 

<f>i{k) = 0(k) = 8^sin 2 (vrA;,). (62) 
e=i 

Given a state (k, i) at t = 0, it jumps at time r to the state (dk', j) with a prob- 
ability i\ t i(dk', j)/0(k), where r is an exponentially distributed random variable 
of mean 0(k) _1 . As before the position process X(t) is defined through 

i r* 

X(t)=X(0) + — / dsVu(K(s)). (63) 
27T i 



6. Homogeneous case: correlations, energy current and 

conductivity 

6.1. Translation invariant measures. We consider a situation where the ini- 
tial measure on phase space is invariant under space translations. For simplicity 
we work in the one-dimensional setting. Using the methods from section [U the 
generalization to d > 2 is straightforward. 

Since energy will be now a.s. infinite, the results of section [3] do not apply. Let 
us denote with <( • )> the expectation value with respect to this initial translation 
invariant measure, and assume that it has the following properties: 

(di) = o, Vy e z, 
(d2) (V(y)V(0)> = 0, £ Z, 

(d3) £^|(VWW*)>|<°°- 

The Wigner distribution is still well defined for every function J £ S(R x T). 
Using translation invariance 




y + y)/2,ky 
J(e(y' + y)/2,k) 



= X>(°)V(*)> / dke-^ kz [-J2Ae(2y + z)/2,ky}, 
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which is finite by condition (d3) and by the fast decay of J. In Fourier space the 
previous expression becomes 

(J,W £ ) 

= ^J2 f d P I ^e 2m(fe - £p/2)y ^(y')>(2/)>e" 27ri(/c+£p/2)?/ ^,fc)* 

and using translation invariance 

(j,yr) = i f dp f dk5(p)W(k)J(p,k)* = - [ dk W(k)J(0, k)*, 
where W(k) is the Fourier transform of the correlation function (■0(O)* , 0(^)), 

W(k) = J2e- 2mkz (mMz))- (64) 

zGZ 

By condition (d3), W(k) is well defined and in Li(T). Moreover, by translation 
invariance, W(k) is a real positive function. 

If we consider the deterministic dynamics only, then W is preserved by the 
dynamics, i.e. dtW(t) = 0. This follows from eq. fHM for 7 = 0. Such property 
is no longer true if the system evolves according to the full dynamics, defined 
trough the generator L = A + e^yS. In order to observe an effective change of the 
covariance, hence of the function W, we have to consider the time scale of order 
e^ 1 . Denoting by W £ (t) = W(t/e), we obtain the following evolution equation, 

d t w £ (k,t) = 1 (cw £ )(k,t) -^[c(y £ + y)](k,t), (65) 

where C is the collision operator defined in ( |22l) . while y £ (k,t) is the Fourier 
transform of the correlation function /■0(O)'0(z)) at the rescaled time t/e, 

y £ (k,t) = J2^ 2nlkz (Ht/e,0)^(t/e,z)). 

As before, we prove that in the limit s — > one obtains a closed equation for 
W £ (t), as stated in the next theorem. 

Theorem 10. Assume that the initial state satisfies the above conditions and 
that W e (k,0) = W (k) is continuous on T. Then, WkeT,te [0,T], 

limW £ (M) = W(M), 

where W(k, t) satisfies the homogeneous Boltzmann equation 

d t W(k,t)= 1 (CW)(k,t), 
W(k,0)=W (k) [ j 

with C defined in |Hj). 

The proof of this theorem will be given in section 16.51 
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6.2. Equilibrium time correlations. We consider the system in equilibrium 

and we denote by ( ■ ) T the average at respect to the equilibrium measure with 
temperature T. This is a translation invariant Gaussian centered measure with 
zero mean, uniquely characterised through its covariance 

W(k)=T, (VW(0)} T = 0, WyeZ. (67) 

( • )„ is a stationary measure for the SDE ( TH1) . 

Consider a function g e i\{Z) antisymmetric, g(z) = —g(—z), and such that 
a; 1 1^ < oo and define the function 



The total time covariance is defined as 

F%t) = Y,(Ht/e)TMO)) T . (68) 

zez 

We want to compute !F E (t) in the kinetic limit e — > 0. 

Consider the centered translation invariant Gaussian measure defined by the 
following covariance, 

>V (T ' r) (*) = rp i "ff. , y (*) = (69) 
1 1 uj(k) + irg{k) 

with r > 0. Observe that W^'^ is a real, continuous function which is positive 
for r small enough. Formally (|69|) corresponds to the perturbed measure 

Z" 1 exp [ - T- X H + rJ2 T M ■ 

We denote by ( • )r TT \ its expectation. 
Lemma 11. For every e > 0, 

£ <*(*/e)T.S(0)> r = lim<$(t/5)> (TiT) . (70) 

The proof of this lemma is given in section [6761 below. 
By direct computation, at the rescaled time t/e, 

< $ (^)> (Tir) = -\^9{z){q z Po- qoPz) {TT) 



uoyk) 



T 



By Theorem [101 W (T ' r) (A;, t/e) -> W(fc,t) for e -> 0, where W(Jfe,t) satisfies 
the homogeneous Boltzmann equation (|66|) with initial condition W(fc, 0) = 
ci;(A;)(T _1 u;(A;) + irg(k))^ 1 . It is easy to verify that for any bounded antisym- 
metric function / on T 

dk (Cf)(k)W(k,t) = - [ dk (f>(k)W{k,t), 
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with = —$(k) = | sin 2 (vr k) [1 + 2 cos 2 (vr k)}, see (El). Then 

/ dk ?^ W (T,T)^ t) = _ % f dk l| W M( M)e -#)t 
Jt v(k) J T u(k) 

and finally, for t > 0, 

lim^(t) = lim-- [ &Mw^(A,0)e-^ 



(71) 



6.3. Energy current time correlation. The Hamiltonian energy current J is 
implicitly defined through the conservation law 

By direct computation J = ^2 z>0 jo, z with 



1 2-1 

)0,z = -^a(z) ^(qz-yP-y - q-yPz-y • 
y=0 



Denoting by ( ■ ) the expectation value with respect to some translation invariant 
centered Gaussian measure with covariance W, it is easy to see that 



[J) = ~\ z a ( z )(lzPo - QoPz) 
o 

-L / dk ^§W(k) 
4:71 J T uj{k) 



2 

z>0 



1 

2^ 



dk u'(k)W(k). 



Let us denote by C e (t) the energy time correlation function on the kinetic time 
scale t/e at temperature T, 

C £ (t) = J2(J(t/e)r x J(0)) T . 
Using the translation invariance of the Gaussian measure ( ■ ) T we have 

C£ w = EEfW E j"CO((?*po - ?oP*)(*/e) 

(q x +z'Px - qxPx+z')(o)) T = Y (^(*/ £ ) r ^(°)) T! 

where 



J = ~7$^ <*(z)(q 2 
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Using the results of the previous subsection we arrive at 

^-w?L^^ m (72) 

dk \u\k)\ 2 e-^ m 



2 



T 

~ 4tt 2 ./ 

for all t G R. 



6.4. Thermal Conductivity. A standard definition of the thermal conductiv- 
ity, kS £ \ is by means of the Green-Kubo formula as 

1 f°° 

k^ = Y2 J dtC £ {et). (73) 

We refer to [21 H] for details. In general, the kinetic limit provides the lowest 
order approximation in e as 

Inserting in (1731 the limit (172]) thus yields 

S «='/*W (74) 
For the pinned case < oo. has been computed in [2j with the result 

«(«o= e -i„(°) +e7 . 

Thus, somewhat unexpectedly, the kinetic theory captures already the main de- 
tails of the conductivity. 

For the unpinned case = oo, hence = oo, for d — 1,2. The Boltzmann 
equation (125]) provides a simple explanation. For small fc, = 27rc|fe|. Thus 
small k phonons travel with speed c. On the other hand, the collision rate vanishes 
as k 2 for small k, see ( 1311) . Thus at small k there are only very few collisions which, 
together with c > 0, is responsible for the divergent conductivity. The positional 
part X(t) of the process consists mostly of very long stretches of uniform motion. 
In fact on a large scale X(t) is governed by a symmetric Levy process of index 
a = 3/2, see [10] for details. 



6.5. Proof of theorem 1101 The proof of the theorem [10] is analogous to the 
proof above. We only have to control that W e (t) and y e (t) are well defined for 
every t e [0, T]. This is stated in the next lemma. 

Lemma 12. Let the conditions (dl-d3) hold. Then W £ {t),y £ {t) G L X (T) for 
every t G [0, T] . 
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Proof. By similar computations as above we find that W E (t), y £ (t) satisfy the 
following evolution equations, 

d t W(k, t) = 7 (CW £ )(A;, t) - ^(c(y £ + y E *))(K t), 

S — 

+ \{c{y e + y £ *m, t) - ^c(w%k, t) + w%-k, t)). 

In particular by Duhamel's formula we can rewrite the second equation as 

pt i A 

37 £ (fc,t)= 7 / ds e- 2luj{k){t - s)/£ (-(3(k)(y £ + y £ *)(k,s) 

Jo 2 

+ \(c(y £ + y £ *))(k, s ) - \c{w £ {k, s ) + w £ (-k, a ))). 

Then we get the following bounds 

|W £ (M)| < W(k) +7C1 / ds [\W £ {k,s)\ + \y £ {k,s 

Jo 

+ 7 c 2 f ds f dk (\W £ (k,s)\ + \y £ (k,s)\) 
Jo Jt 

\y £ (k,t)\<+ 1C3 f ds [\w £ (k,s)\ + \w £ (-k,s)\ + \y £ (k, s )\] 

Jo 

+ 7 C4 f ds f dk (\W £ (k lS )\ + \y £ (k,s)\) 
Jo Jt 

and finally 

dk [\w £ (k,t)\ + \y £ (k,t)\] < [ dk\v(k)+ 

Jt 

7 c 5 f ds I dk (\W £ (k,s)\ + \y £ (k,s)\). 
Jo Jt 

By Gronwall's lemma 



dk [\W £ {k,t)\ + \y £ {k,t)\] < e 7C5 ' / dk W(k). 
t Jt 



□ 



6.6. Proof of (1701) . We define the generator of the speeded up process as 

L £ = e~ l L = e~ x A + S. 

Let us denote by ( • ) T the average with respect to the Gaussian measure with 
zero mean and covariance (1671) , and by ( • ) . the Gaussian measure with zero 
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mean and covariance (|69|) . We consider the Laplace transforms of ($(t/e)), T 
and E«6z(*(*/e)T,$(0)> T , 



/CO 

* (t/e)) (TjT) = <(A - L,)" 1 *)^ = < MA ) (Tjr) , 

/CO 
dt e"* (Ht/e)r x m) T = £ <[(A - ^) _1 *M> 



(75) 



' T' 



where «a = Ezez <?o£>,z with fx the solution of the equation 

A f x (z) - ^A(4/ A («) + / A (z + 1) + / A (* - 1)) = </(*). 
6 

Observe that u\ does not depend on e, because L £ u A = e~ x V F + 7 Sua with F 
some non-local function. Since /a is antisymmetric, and by translation invariance 
of the measure, the gradient term does not contribute. We have 

( Ux )t,t = E /a(«)(?qP*>t,t = ~~\ E f^K^Po - QoPz)t,t 

■it* m 1 



T A + 7 0(A;)T- 1 w(fc) +ir£(fc) 



with 0(/c) = —(3(k). For every positive A, the right hand side is finite for every 
r. In particular 

1/.. v l«(*)l J 



In the same way 



t w(*) a (A + 7#A0)" 



E <^ r.*> r = T E E /A(y)p(«)<9b?»-.)> 5 

™rf/JW^ _ ^2 / » lff(fc)| ; 

where Parseval's identity is used. □ 



T [ dkUky^%W^\k)=T* [ dk 

Jj w(ky jj 



7. Appendix 

Proposition 13. Under the assumption (b3), 

\(J,W)\<Cj (76) 
for every test function J G S(M. d x T d ) and for every e > 0. 
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Proof. For any test function J G S(R d x T d ) we denote by || J||„ i00 the following 
norm, 

d d 

\\J\\n,oc= sup sup | V... V<9 fcii ...<9 fcin J M -(r,k)|. 

i,je{i,...,d} reRd keT< i h=1 

Let us define 



Ji d 



for every y, y' e i,j = l,...,d. Integrating by part in k for (d + 1) times, we 
get 

\mJ](y,y',i,i)\ < (2 ^ + .| fa; -, l) + .. 1 +(;/ ;- ;/J )| J+1 ll J ll^- 

where t/i denotes the i-th component of the vector y. We denote by 

(y'-y) d+1 = \(y , 1 -yi) + --- + (y'd-y d )\ d+1 - 

By Schwarz inequality 
\(J,W £ )\ 

d 

<(e/2) d EKl^(y)l 2 >J 1/2 Kl^(y')l 2 )J 1/2 |W £ [J](y,y / ^,j)l 

y,y'£Z d i>i=l 



<(s/2y(\m 2 ) E J2 7 ^mCo\\J\U + i,oo< 



cy\\d+i,c 

d 



where in the last inequality we used £ d (||?/>|| 2 ) < K with K positive. □ 



Lemma 14. Let Assumption[J\hold with a(0) = (unpinned case). The following 
assertions hold. 

(i) There are constants C\, C2, C3 such that Vk £ T d 

|Vd(k)| < d|k|, C 2 |k| < w(k) < C 3 |k|. (77) 

In addition, ||Vu;||oo < °o. 

(ii) For all k e T d and p G M d ; i/iere is a positive constant C such that 

e _1 |w(k + ep/2) -w(k-ep/2)| < C|p| (78) 
for every e > 0. 

(iii) For p e M d , i/iere zs a positive constant C4 stzc/i t/iat for all k e T d , 
|k| > e|p| ; with e > one /ias 

Inl 2 

\e- l [uj(k + ep/2) -u(k-ep/2)}-p- Vw(k)| < eC 4 ^- . (79) 
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Proof. The first inequality of item (i) follows from a Taylor expansion of d 
around zero, using the fact that Vd(0) = and H-D^Hoo < oo. Using the same 
argument we have uj(k) < C 3 |k| , since u;(k) = d(k) 1//2 and d(k) < C|k| 2 , with 
C > 0. Let us denote with Aq the Hessian of at at k = 0. By assumption 
(a4), k • Aok > c|k| 2 . Moreover, there is a 5 > such that for every |k| < 5, 
(|d(k) - |k ■ A k\) 2 < ±k • A k. Then w(k) = ((d(k) - |k • A k) + |k ■ Aok) 1 ^ > 
||k| 2 if |k| < 5. For |k| > 5, uj is strictly positive and there is a constant c' such 
that u;(k) > c'|k| if |k| > 8, which proves the second inequality. To prove the last 
one it is enough to observe that for all k ^ 0, 

_ 1 |Vd(k)| Id 

Hence ||a>||oo < °°- 

Let us prove item (iii). We observe that the function on the left hand side 
of (1791) is zero if |p| = and we have to consider only the case |p| > 0. Since 
we are assuming |k| > e|p|, it follows that we need to discuss the case |k| > 0. 
Observe that for any s G (0, e], |k±~sp| > ~|k| > if |k| > e|p| and the function 
u(k + \ey>) — u(k — |ep) is C°° in this range. In particular 

e _1 u;(k + ep/2) - u{k - ep/2) - p • Vcj(k) 

= e (| • V) 2 w(k + sp/2) -e(f.v) 2 .(k- 5p/2) 
with s,s£ (0, e), where, denoting k ± = k ± |sp, 

(p • V) 2 u(k ± ) = (p ■ V) 2 d(k ± ) - (p • Vd(k ± )) 2 

and, using item (i), 

|(p-V)Mk ± )| <4C 4 ^. 

|k| 

This prove item (iii). 

Item (ii) follows from (iii) if |k| > e|p|. If |k| < e|p| we use the bound 

£~V( k + ep/2) - u{k - ep/2)| < C 3 £-\\k + ep/2| + |k - ep/2|) < C|p|. 
7.1. Proof of ( 1381) . First of all we observe that 

+ [Y y ^ y 'y}[Y y ^(y)} + [Yy^MyTWy-My)}, 



10) 



where 



Yy+My) = 77|(Py+i ~ Pv+i)> Y vi>{y) = ~ Pv-*' 

Y y -ilp(y) = -j=(Vy-1-Vy-\)- 



51) 
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Then f lHUj) is equal to 

(Py+l - Py+2) [{Py+l ~ Py+2)Sy',y + (Py+2 ~ p y )6 y > >y+1 
+ (Py - p y+1 )5 y/ty+2 ] + (Py+l ~ Py-l) [(p y ~ Py+l)5y> ,y-l 

+ (Py + l - Py^l)5y l} y + (j)^ ~ Py) Sy> >y+1 ] + (p y ^ 2 ~ Py-l) 

X [(Py-l - Py)$y>,y-2 + (Py ~ p y - 2 )Sy', y -i + (p y - 2 - Py-l) <V ,y] 

= (Py+l - Py+2)(Py - P y+ l)5y>, y+2 + (Py-2 ~ Py-l)(Py-l ~ Py)$y>,y-2 

+ [(Py+l - Py+2)(Py+2 ~ Py) + (Py+l ~ Py-l)(Py-l ~ Py)} $y>,y+l 

+ [(Py+l ~ Py-l)(Py ~ Py+l) + (Py-2 ~ Py-l)(Py ~ Py-2)] V.v-l 

+ [iPy+1 -Py+2) 2 + (Py+l -Py-lf + (Py-2 - Py-l) 2 ] <V,2/ 
2 

= 2 A y>r S y ' )y+r 

r=-2 

and we can write 

\{e/2) ([YMyr][YMy)]) e J(e(y' + y)/2,y- y r 

y,y',z& 

= l(e/2)J2i2 A v^(y + r/2),-ry. 

° yeZr=-2 

Expanding the A y>r and identifying terms that differ by translations, we arrive at 
1(8/2) J2 {[YMyyWMy)]) £ j{e(y' + y)/2,y-yy 

y,y',zeZ 

= f E [< 2 »+l-»+2-^> e (^' 2 )* + ^N'- 2 )*) 

y& 

+ (2p y p y+2 - 2p 2 y ) £ (J(ey, 1)* + J(ey, -1)*) 

+ < - 4p y p y+ i - 2 PyPy+2 + 6p 2 y ) £ J{ey, 0)*] + 0{e), 

where we have used the smoothness of J in x G M. We can rewrite the last 
expression as 

2 

( e / 2 )^ Yl ®(z,u)(p y p y+z ) £ J(ey,u)* + O(e), 

y£X z,u=—2 

where a(z, u) = a(z, —u) = a(u, z) = a(—u, z) and is given by 

a(0, 0) = 1, a(0, 1) = -1/3, a(0, 2) = -1/6 

a(l,l) = 0, a(l,2) = l/6 (82) 
a(2,2) = -1/12. 
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